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Abstract 

A Langevin equation with a special type of additive random source is consid- 
ered. This random force presents a fractional order derivative of white noise, 
and leads to a power- law time behavior of the mean square displacement of a 
particle, with the power exponent being noninteger. More general equation 
containing fractional time differential operators instead of usual ones is also 
proposed to describe anomalous diffusion processes. Such equation can be re- 
garded as corresponding to systems with incomplete Hamiltonian chaos, and 
depending on the type of the relationship between the speed and coordinate 
of a particle yields either usual or fractional long-time behavior of diffusion. 
Correlations with the fractional Fokker-Planck equation are analyzed. Pos- 
sible applications of the proposed equation beside anomalous diffusion itself 
are discussed. 
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1 Introduction 



In recent years, growing attention has been focused on the processes that 
take place in random disordered media, and in dynamic systems demon- 
strating chaotic behavior. A special place here is taken by the systems with 
incomplete Hamiltonian chaos, which trajectories in the phase spacecan be 
portrayed as a set of "islands around islands" with self-similar (fractal) struc- 
ture. In such systems the lifetime of any state of regular motion is also ran- 
dom. Incomplete chaos results in a few interesting phenomena that, among 
others, include anomalous properties of transport processes. That is asso- 
ciated with the fact that the islands of stability in this case become to act 
as a system of traps with a certain given distribution of the trapping time. 
For example, well known is the phenomenon of anomalous diffusion that is 
characterized by the time function of the mean squared displacement of a 
diffusing particle, which is described not by the Einstein's law but by an 
power function with fractional exponent [0, |^ . 



In most cases, such behavior is considered to be connected with self-similar 
properties of the diffusion medium. As this takes place, the Fokker-Planck 
equation describing such diffusion process was shown to involve the integrod- 
ifferential Riemann-Liouville operators I'^ and of fractional order (see 
Appendix A). Possible interpretation of the physical meaning of these oper- 
ators was offered in |^, and assumes that a system described by equations 
with fractional derivatives or integrals possesses a "selective" memory that 
acts only in the points within a set of dimensionality u, and is in accordance 
with the ideas about a certain self-similar (say. Levy-type or fractal) distri- 
bution of traps and waiting times Fractional Fokker-Planck equations 
were studied in 0-0, and proved to be quite useful to model anomalous 
diffusion processes. The most general form of such equations was introduced 
in H, H and in a simple case looks 



where n is the concentration of diffusing substance and the fractional dif- 
ferential operators on the right-hand side and on the left-hand side are ei- 
ther fractional Riemann-Liouville derivatives or so called local fractal 
derivatives 0. 
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Another method of the anomalous transport properties description refers 
to their microscopic motion and involves introducing fractional derivatives 
into a stochastic process that characterizes the Brownian motion. As the 
starting point to be changed the Wiener stochastic process that defines the 
displacement of a particle with time is taken 



X 



t 

Xo = j F{t')dt' 



where F{t) is the the Gaussian white noise (the brackets here mean averaging 
over the possible realizations of the process), 

and the fractal Brownian motion is obtained by substituting the Riemann 



integral for a fractional one |T0|-p!3| 





The fractional integral of the white noise was named the fractional noise. 
This noise is Gaussian but nonstationary. Its correlation function has the 
form 

(x(ti)a;(t2)) = 2Fi{l, 1 - a; 1 + a; t</t>) (3) 

a [1 [a)\ 

= min{ti,t2), t> = maa;(ti, ^2) 
(x(t)x(t)) oc t'"-^ 

In principle, this noise can be made stationary, but only provided that —00 
is taken in Eq.(^) as the lower limit of integration instead of 0. However, in 
this case one has to somewhat modify the power kernel in order to provide 
the convergence of the integral. We will briefly discuss this below, at the end 
of the Section 2. 

Fractal Brownian motion is of great interest not only from mathematics 
and theory of stochastic processes points of view but also in terms of physical 
applications. It can be used in describing polymer chains, electric transport 
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in disordered semiconductors, diffusion on comb-lilce structures, and so on. 
However, it is well known that in many cases the most convenient way of 
describing the Brownian diffusion of particles is not the Wiener process but 



rather the Ornstein-Uhlenbeck process (or Langevin method p|, 0, |T5|) that 
is based on the solution of stochastic differential equation 

jv = -IV + Fit) (4) 

where v is the velocity of a Brownian particle, 7 means the factor of liq- 
uid friction, and F{t) is the random source characterizing the properties of 
medium where diffusion occurs. Displacement of a particle and its path are 
determined not directly but through the integration of the instant velocities. 
Apart from the problem of Brownian motion itself, this method is widely ap- 
plied to describe various systems subjected to external noise. This brings up 
the question of whether the Langevin equation can be written for anomalous 
diffusion as well, and if so, what will be the structure of the corresponding 
random source in it. In principle, if the Fokker-Planck equation is known 
then the Langevin equation can be derived from it. However, this method 
involves certain difficulties, and gets even more complicated because of the 
intricate structure of fractional integrodifferential operators. So, we shall 
postulate the form of such equation in further sections of our paper. Actu- 
ally, we shall consider three forms of such equation that have in common the 
method of introducing the memory by means fractional derivatives yet differ 
in the properties of medium and in the ways of setting up the problem. 



2 Langevin Equation with Fractional Deriva- 
tives 

a. As it has been mentioned above, introducing the fractional differential 
operators into the Fokker-Planck equation makes it possible to describe the 
anomalous transport process quite correctly. Therefore, let us consider the 
following equation that differs from the usual Langevin equation by replacing 
the first derivative with respect to time by the fractional derivative of order 



= -^v + Fit) (5) 
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Applying the fractional integral operator to both left-hand and right-hand 
sides of this equation, one has 

v-vo=oIti-lv + F{t)) (6) 

Expressing the fractional integral in the explicit form, we rewrite it as follows 

v = vo + A(t) ^ / — dx', A(t) = —— / — dx' 

^ ' r(z/) J {x- x')^-^ r(z/) J {x- x'y-'' 



(7) 

This equation can be easily solved by standard techniques for Volterra inte- 
gral equations, and its solution has the form 

t 

V = voE^,,i-^n + j F{t'){t - tr^'E,,4-^it - t'YW (8) 



where Ea^jsiz) is so called Mittag-Leffleur function | T6[ 

oo 



If u = 1 then Eq.(H) reduces to the solution of the usual Langevin equation 




V =Vo 



If F{t) is taken, as usual, to be a Gaussian (5-correlated source with zero 
mean, then the velocity correlation function has the form 

{v{tMt2)) = vlE^A—ft'i)E^A-^r^) 

min{t\ ,t2) 


The integral from the right-hand side can be taken in the explicit form only 
if we introduce t> and t< from Eq.(^, which is rather inconvenient. How- 
ever, we can obtain several important results without having to calculate 
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the integral directly. It is only essential that the integral is a symmetrical 
function of its arguments ti and ^2- Further, let us find the corresponding 
mean squared displacement of the particle, which motion is described by the 
equation Eq.(|^). It follows the expression 




{v{ti)v{t2))dtidt2 







{{Axf) 



+ 



t - 2tE. 



2mK 



-in 



1 2 



dt 



(10) 



which, as shown in the Appendix B, for correlation function of Eq.(^) can be 
reduced to 



At large values of the argument Ei^—'-^V) oc E2^u{—lt'^) oc l/'yt'^. Therefore, 



the integral in Eq.(|TT|) necessarily converges at t oo for 1/2 < i/, and the 
leading term in the asymptotics is 



((Ax)2) 



(12) 



just as it is the case for the classical Langevin equation with the first deriva- 
tive with respect to time. Even when v <l/2 and the integral diverges, it 
increases more slowly then the first term and the linear asymptotics remains. 
Leaving the discussion of this point for later, we shall try to modify somehow 
the initial equation. 

b. First, it must be noted that in Eq.(|ll|) we assumed that, as usually, 
the velocity v was defined as the first derivative of the coordinate with respect 
to time and, therefore, x = j^v{t)dt. Actually, it makes sense to consider a 
more general relationship 



X 



v{t') 



mj {t-t'Y- 





-dt' 



(13) 



that corresponds to the most complete possible description of the system 
memory by means of the fractional integrals - the particle's displacement is 
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defined by velocity only in the points within a time interval of dimension 
V. Though in this case the physical meaning of the corresponding inverse 
definition of velocity as a fractional derivative of the coordinate with respect 
to time needs to be explained. In order to clear this point, let us recollect 
that microscopic motion of a diffusing particle represents a twisted and ev- 
erywhere nondifferentiable curve. For such curve, however, one can often 



find a derivative of fractional order that, in fact, is a derivative of the 



trajectory, averaged with a power weight, and the observed motion of the 
particle is the thus averaged motion. In terms of memory, it means that 
for fractal paths (which are even the paths of classical Brownian particle) 
some of the instant velocities and displacements do not contribute into the 
resulting macroscopic motion. As this occurs, the behavior of the solution 
changes, and diffusion becomes anomalous. In this case instead of Eg . (|TTD 
we obtain 



72 Z^^T{v + vk)V{vl) T{vk + vl + 2v) ^ ^ 

At large times the first term tends to ^^"^ series in the second term 

converges (see Appendix C). Therefore we finally have 

{{/\xf)^N\e^-\ 1<N<-^ (15) 

that coincides with the asymptotics obtained in p|, |^ for the generalized 
equation of fractal diffusion Eq.([^) if 7 = 1, and with the asymptotics of 
fractional Brownian motion [|10[-|1^. The larger u the faster the particle 
moves. More accurate estimations of factor A^ can be obtained either follow- 
ing the method similar to that stated in the Appendices B and C (which in 
this case would be quite a challenge) or numerically. 

c. Let us now investigate separately the influence of memory and fractal 
behavior on the regular and random components of the force acting on a 
particle. It turns out that the same fractal asymptotic behavior can be 
demonstrated by taking into account only the memory for the random force 
component in the initial Langevin equation, which has a more clear physical 
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meaning (cf. Eq.(|^)): 



f = f — 7 / v{t')dt' + 



m 
it - t') 



-dt' 



(16) 



In fact, it means that in Eq.(^ the random source is not 5-correlated but 
represents the fractional derivative of the 5-correlated process g{t), which can 
be taken in the sense of the generalized functions |jl8| 



F{t) = oDl^git) {g{ti)g{h)) = q6{t, - h) 



(17) 



Similar equation was actually considered in ||20[. It was shown numerically 
there that when in the Langevin equation Gaussian white noise is replaced 
by the fractional Gaussian noise JTO], |11[] it can yield the spectra for homoge- 
neous Eulerian and Lagrangian turbulence. However, equation Eq. (p!6D can 
be easily solved analytically, and the solution reads 



z 

Voe-^' + J F{t'){t-tr-'E,,^[-^{t-t')]dt' 



(18) 



The corresponding expression for correlation is 



min{ti ,t2) 



+ q 



dt 



(19) 



Calculation of the mean squared displacement is similar to that in the Ap- 
pendix B and results in 



2^ ..(l-e-T 



+ 4^ 



s 



-7) 



k+l 



k+l+2iy+l 



k,l=0 



+ V + k)V{u + l){k + l + 2p){k + l + 2p+l) 



(20) 



By differentiating the sum two times, we obtain 



dt^ ~ dt 



fEi+,^,{-^t) - 1 
8 



from which it follows that when t ^ oo it is the second term that defines the 
asymptotic behavior of the displacement 



(2z/-l)[r(z/ + l)]'7' 



((A^)') « ^t::^ 77724^^'^-' (21) 



which agrees with Eq.(|T5D up to the constant multiplier. 

d. Finally, let us consider the third case, that is when the memory is 
taken into account only for the friction force 

t t 



where F{t) is again the Gaussian white noise (cf. Eqs.(|^), (]T6D). It means 
that now it is the dissipative force that is proportional to the fractional 
derivative of velocity fdiss = oDl^^vit). Its solution is 

t 

V = voE,,,{-^n + J F{t')E,^, [-7(t - t'Y] dt' (23) 



The contribution from the second term in the mean squared displacement 
results in the power function of time in the following form 

{{Axf)^M\t'-'\ 1<M<—^ (24) 

Y 1 (1 + ^) 

It should be noted here that, unlike the previous two cases, here no limitations 
are imposed on the value of z/, and depending on whether it is greater or 
less than 1, we can obtain either subdiffusion with the power exponent less 
than 1, or super diffusive behavior with the power exponent greater than 
1, observed, for instance, for a phase for wave propagation in nonlinear or 
random media. The greater is u, the slower the particle moves (it was the 
other way around in the previous cases), because at z/ > 1 the corresponding 
dissipative force is no longer a fractional derivative of the velocity but a 
fractional integral, and, hence, its influence increases. At z/ > 3/2, the 
mean squared displacement Eq. (|2^) tends to zero with time increasing, which 
corresponds to a stop in the particle motion since the terms neglected in 



Eq.(gJ), as well as in Eqs.(14)- ([r5|) , are either vanishingly small or constant. 
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and means that particle energy dissipates faster than it is pumped. On the 
other hand, hmitations 1/2 < u < \ were used in the previous cases only to 
simplify the estimations for the sums which were shown to be convergent and 
therefore in principle it must not affect the main results and they also can 
give sub- and super diffusion behavior. But because of the fact that prehistory 
influences only on the dissipative force, now diffusion is impeded when we 
increase the value of v. 

Let us now analyze the results obtained. If we substitute the first deriva- 
tive in the Langevin equation for the fractional one but use the same relation- 
ship for velocity and coordinate, then the solution possesses the same linear 
asymptotic behavior as the initial solution does. This means that taking 
into account the memory for the friction force and random force at the same 
moments of time does not affect the particle's motion at large times but only 
in the beginning of the motion. In a sense, this is quite reasonable, since 
a random source provides a particle with an additional energy, and friction 
results in its dissipation. At larger times, when the system arrives at the 
stable state, these two processes compensate each other, but only provided 
that they have the same duration. 

Using the relationship Eq.(|T3|) means that the particle's motion equation 
written in the Newtonian form reads 

and is the most general in the sense of memory calculation - every deriva- 
tive with respect to time becomes fractional. Then anomalous properties of 
the diffusion arise from fractional integration of velocity, that is, from the 
assumption only part of instant velocities contribute to the final path. 

The anomalous diffusion is also the result of the Langevin equation with a 
source that is a fractional derivative of the white noise. The autocorrelation 
function of such noise, just like that in Eq.(^, has a power behavior yet 
negative values of the exponent and therefore appears to be a generalization 
of the flicker-noise (see 0]). It should be noted that a stationary noise can 
be obtained only by using the Liouville fractional derivative with the infinite 
lower limit. However, when used directly in the equations Eqs.(0), (|TBp, 



p^), these derivatives would lead to the divergences because of the power 
integrands. Besides, this is not quite reasonable, since in real situations there 
is always a moment at which the motion starts, and assuming this moment as 
infinitely distant is possible only in a system where characteristic relaxation 



10 



scales exist. Fractional Brownian motion, however, can not be assigned to 
such processes - its memory is described by a power function and leads to 
the long-range correlations having no time scale of their own. Nevertheless, 
the nonstationary of the noise should not lead to a misunderstanding. As we 
shall illustrate below, such process can be regarded as stationary only in a 
wide sense (the corresponding transition probabilities depend uniquely upon 
the duration of the transition At) only for small At, that is for At << t. 

The absence of stationarity is an important and rather obvious prop- 
erty of anomalous diffusion, although it is not always given proper consid- 
eration. Here it is relevant to note the following. Normally, the Brown- 
ian motion can be described with two different stochastic processes. The 
Ornstein-Uhlenbeck process is strictly stationary but does not have indepen- 
dent increments. Moreover, its increments are not even uncorrelated. The 
Wiener process, which is the integrated Ornstein-Uhlenbeck process in the 
limit of intense friction and noise, has stationary independent increments, 
but is neither strictly stationary nor a wide-sense stationary. The process 
under investigation with a fractional derivative of the white noise seems to 
be an intermediate process that does not have stationary increments but is 
asymptotically stationary in a wide sense. The fractal nature of motion leads 
to the fact that even this kind of stationarity is observed only within short 
periods of time, which means that the process becomes quasistationary. This 
conclusion is in a good agreement with the notion that anomalous diffusion 
is just an intermediate asymptotic behavior for systems of certain types. 

Finally, if the prehistory affects only the dissipative force acting on a 
particle, its behavior also becomes anomalous. 



3 Probability distributions 

Let us now study the probability distributions that arise from the stochastic 
processes described above. For the sake of convenience, let us consider only 
the process with the noise that is a fractional derivative of the white noise 



Eq.(|T^). The other cases are considered in perfect analogy to this one and 
give similar results. 

The distribution of the particle coordinate will be Gaussian. It follows 
from the fact that the equations are linear both with respect to the external 
additive noise and the stochastic variable. In order to determine the transi- 
tion probability P{x, t + At; y, t) let us consider transition moments at large 
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times 



t+At t+At 



Mn = \x - y\ 



/ / \ / \\ 7 7 ) n = 2k + I 

{V{ti)...v{tn))dti...dtn = < (2fc)! ^ _ 

(26) 



where 



t+At t+At 



{v{ti)v{t2))dtidt2 



t t 

t+At t+At min{ti,t2) 



dti I dt 



This imphes that the probabihty of transition obeys the Gaussian distribution 
with dispersion I. In order to identify time dependence of the dispersion, let 
us rewrite the internal integral in the form /q*^^* = /J + jl^^^ ■ In this case 
Eq.(^) will have two terms, the first of which allows changing the order of 
integration since the limits of integration are no longer dependent on the 
intrinsic variables: 



t r *+^* 



dt' 



^ ..1 [-7(^1-0] 



dti 



At At << t we then get 
q At^ 



A 



72 t3- 



2u' 



2u[T{u)\- 



<A< 



(28) 



The second terms that appeares in Eq.(p7D after the substitution ti — ti — 
t, t2 ^ t2 — t, t' f — t reduces to the mean squared displacement given by 
Eq.(^) during time At 



(2t.-i)[r(z.+ i)]^^ 
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which is much greater than I2 if At << t. Hence, 



'72(21/- i)[r(t/ + 1 



X exp 



2Tiq Af^^-^ 
7^(2z.-l)[r(z. + l)]' 



X 



y\ 



2q 



At 



2u-l 



(29) 



This expression also gives a "suitable" result for the coordinate distribution 
W{x, t) of the wandering particle that agrees with 
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W{x,t) = P{x,t;0,0) 
B 



1 

V27r5t2^-i 
1 q 



X 



exp 



(30) 



2 0,2 



(2z.-l)[r(z. + l)]^7 

and follows from the fact that if we take ?/ = 0, to = as the initial point of 
transition, then the I2 = holds disregarding the value of At. Therefore, at 
large times the described process becomes quasistationary in a wide sense, 
that is, the probability of transition during a small time interval At depends 
on the value of At only. However, if the condition At << t is not fulfilled, 
then I2 can no longer be given with the simple expression like Eq.(^) and 
thus can not be neglected as compared to Ji. It should also be noted that such 



process on any time scales is not Markovian since the equations Eqs.(29)-(|30D 
do not comply with the Chapmen- Kolmogorov equation 

P{x, t + At) = J P(x, t + At; y, t)P{y, t)dy 

which can be verified by direct substitution. Indeed, if we introduce in 
Eq.(^) new time variable as r = t^^~^, then the Markovian Gaussian process 
will be that with the transition probability distribution 



P(a;, r+Ar; r) 



and at small At 



V27rL>Ar2'^-i 



exp I 



\x - y\ 



2DAt 



2u 



At ^ {2iy - l)t"'-'At 



2u-l ,2u-l 



- Ar = (t+At)2'^-^-t 
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It is worth to note that in this case the property of wide- sense stationarity 
is lost for all time scale. 

Non-Markovity is a consequence of the memory of the past process exis- 
tence. If for a Markovian process the future is uniquely determined by the 
present, then in the initial equation Eq. ([T6|) the behavior of a particle in the 
next moment of time is, generally, dependent on the whole previous history 
starting from the very beginning of the motion, and the transition probability 
P{x, t + At; t) does not depend on the time t only when At — >• 0. In this 
connection the question of transition from Eq. (|16]) to the Fokker-Planck-type 
equation still remains to be solved. Indeed, the proposed in 0, |[ method 
is based on the fractional Taylor series expansion of the transition probabil- 
ities in the Chapmen-Kolmogorov equation. For example, according to 0, 
the fractional Fokker-Planck equation for the transition probability given by 
Eq.(^) takes the form 

Q2U-1 Q2 

^^W{x,t) = D—W{x,t) (31) 
Solution of this equation was obtained in [|, ^ and shown to be 

^'^^''^ - Vm^ (1,1) ) ^^^^ 



where HI'^ is the so-called Fox function ||T7| , p^ . However, according to 
Eq.(ig) 

W(x,t) = , ^ exp ( = H''^ ( —^\ , " , 1 ^ Wi(x,t) 

(33) 

Therefore, the solution Eq.(|55D does not satisfy Eq.(^T]) though the time and 
the coordinate of a particle appears in its solution in the "correct" combi- 
nation x"^ /t^^"^. However, as the stochastic process under consideration is 
not globally stationary, this fact seems to influence the structure of the cor- 
responding Fokker-Planck equation. In particular, such equation may also 
include the terms that characterize the sources and drains of probability, and 
therefore its structure will no longer be as simple as in Eq.(^). Probably 
this problem may be solved by studying the waiting time probabilities and 
then writing down the equations similar to those used in 0. 
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The average kinetic energy of a particle can be obtained as 



(E) = {[vm (34) 

Substituting here the expression for the velocity correlation Eq. (|TU|) and tak- 
ing into account that at large times, when the motion becomes stable, the 
particle "forgets" the initial velocity and we can neglect the first term in it 
(see also Eq.(^)), yields 

t 

{E) = ^lr'^'^E.,,{-jr)fdx (35) 



Since when t — oo E^^i{—'-yT) oc I/7T and in the vicinity of zero the 
Mittag-Leffleur function is bounded, for 1/2 < v the integral in Eq.(P3D 
converges. Therefore (see [|14|, the velocity of the moving particle will obey 
the Maxwellian distribution. 

Probability distribution for the processes described in the paragraphs 2b 
(Eqs.(^ and ([T3|)) and 2d (Eq. (|22|) ) are obtained in the same way. Veloc- 
ity distribution obeys the Maxwellian law, and the coordinate distribution 
W{x^ t) and transition probability P(a;, t + At; y, t) at large times t and small 
At follow the Gaussian law with the dispersion determined by the mean 
squared displacement over times t and At, respectively. 



4 Conclusion 

To conclude, we would like to note the following. The equations and models 
proposed in this paper are quite general and, represent a way to introduce 
and describe a certain class of Gaussian non-Markovian stochastic processes. 
Processes of this type are featured with the "selective" memory acting only 
in the moments of time distributed over a Cantor-type fractal set and taken 
into account by means of fractional derivatives, and show promise in de- 
scribing stochastic processes in fractal media and systems with incomplete 
Hamiltonian chaos. One of such processes is anomalous diffusion, observed 
in a wide variety of systems. The method developed in this paper is an al- 
ternative to that utilizing the fractional Fokker-Planck equations, although 
we could not trace fully the relation between these two approaches and it 
still remains to be explored. Nevertheless, the two approaches give similar 
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results, and choosing between them should probably be determined by a type 
of problem to be solved and by practical convenience of calculations. Earlier 
the fractal Brownian motion was represented only as a Wiener stochastic 
process of fractional order (the integral of order u from the white noise or, 
which is the same, the integral of the first order of the noise which is a frac- 
tional derivative of order z/ — 1 of the white noise). As shown in the present 
paper, the Ornstein-Uhlenbeck process with fractional noise leads to similar 
results and allows to obtain more general probability distributions in an eas- 
ier way, as compared to the path integrals [0. The equations proposed can 
also be easily generalized by adding linear and nonlinear terms, and be used 
for calculation of various statistic characteristics of real stochastic dynamic 
systems. 
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A Fractional dervatives 



There are about two dozens of different definitions for fractional derivatives 
that are in one way or another adapted to various features of classes of 
functions for which they are defined. The most comprehensive description 
of the mathematical aspects of the issue is given in the monograph In 
physics the Riemann-Liouville fractional derivative is the most commonly 
used. Its definition goes back to the well-known Cauchy formula for multiple 
integrals 



J- fit) = J dt...j dtfit) = ^^^^ J 



{t - V) 



Substituting the factorial for the Euler gamma-function, we can generalize 
the formula by introducing the fractional exponent a as follows 

t 

a 

This expression is referred to as the Riemann-Liouville fractional integral. 
The fractional derivative is then defined as an ordinary derivative of the 
integral of fractional order 

1 d"" f fit' 



Tin - a) dV' 



[ dt' n-l<a<n (A.2) 



Fractional derivative can also be treated in the form of convolution with 
power function, and in this sense the definitions in Eqs. (|A.l|) -( |A.2|) can be 
easily transferred onto generalized functions. 
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B Calculation of the mean squared displace- 
ment 

The mean squared displacement of a particle during time t is defined as 



t t t 
{{Axf) = [ [ {v{h)v{t2))dhdh = v^J [ ^i,.(-7t' ndt'] ' + (B.l) 







where the integral in square brackets in the first term equals to tE2^u{~lt''), 
and the second terms takes the form (we expand the Mittag-Lieffleur function 
in a series) 

t t min(ti,t2) 



k 



''-^ 



Since the triple integral here does not change its form when interchanging 
ti ^2, then, assuming for certainty ti < t2, we can write it as 

t t2 h 

4; = 2 y dh j dt, j dt'ih - tr'-\t2 - tY~' 



This integral is now easy to calculate using the rules of fractional integration 



and the properties of Gauss generalized hypergeometric functions 2F1 1T7 

t t2 

hi 



'^jdt2j dhT{uk\^^{t2-hy'-^ 

dt2tf^'' [ C:^Y' 2F,{1, + Vk- 

J ^t2^ 12 12 





t t2 

2T{vk) 



T{vk + l) 




t 



dt2tf+''' = ——. ) (B.3) 



iyk{uk + iyl)J ^ uk{iyk + iyl){uk + + 1' 
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Substituting in Eq.(|TT|) yields 

oo 



-7t 



v\k+l 



r(l + uk)T{ul){uk + ul){uk + z// + 1) 



(B.4) 



In order to evaluate the asymptotic behavior of this series at t — * oo, let us 
note that 



df^ ^ rri + vk\ ^ 



-It 



v\l-\ 



d 
It 



k,l=l 

oo 



r(i 

k=i ^ 



uk) 



d 
It 



iE,A-in - 1 



(B.5) 



Integrated twice, and taken into account that I{t = 0) = ■;^|t=0) for the mean 
squared displacement it gives 



+ 



t - 2tE2,u{-lf] 



dt 



(B.6) 



C The sum evaluation 

We can find the asymptotic behavior of the double sum in Eq. ([1^ (we desig- 
nate it as 5*1) in the following way. First, let us find the upper boundary for 
this sum. To do this we substitute the second cofactor (z/ > 1/2, /c, / > 1) 

T{vk + + v -I) T{uk + + 



T{uk + + 2v) {uk + vl + u- l)T{uk + ul + 2u) 

uk + ul + u — 1 vl ^ ^ 

hence 

(C.2) 
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Therefore, S*! < 1 /T{y) when t ^ oo. In order to get the lower approximation 
for the sum (for l/2<z/<l)wc shall multiply it by t and differentiate twice. 
Reducing the argument in the gamma-function in the numerator hy u < 1 
and increasing the denominator by (2 — 2u) < 1 we then get 



{Sit) > 



-7) 



k+l 



T(iyk + ul - 



£ + yk)T{ul) T{vk + i/Z + 2) 

1 d i-^fr^ i2 Id 



2dt 



^r(l + i/A;)J 



2dt 



Ei,,{-^n - 1 (C.3) 



Since this inequality holds for every t, we can integrate the left-hand and 
right-hand sides, which yields 5*1 > 1/2 for t — > oo. 
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